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Abstract 



We construct a replica field theory for a random matrix model with logarithmic confinement [K. A. Muttalib et.al., 
Phys. Rev. Lett. 71, 471 (1993)]. The corresponding replica partition function is calculated exactly for any size 
of matrix N. We make a color-flavor transformation of the original model and find corresponding Toda lattice 
equations for the replica partition function in both formulations. The replica partition function in the flavor space 
i^i - is defined by generalized Itzikson-Zuber (IZ) integral over homogeneous factor space of pseudo-unitary supergroups 
SU(n | M, M)/SU(n | M — N, M) (Stiefel manifold) with M — > oo, which is evaluated and represented in a compact 
CNI i form. 
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I. INTRODUCTION 



Recently a considerable interest appeared towards analyzing of integrable structure behind Gaussian random matrix 
theories (RMT) 1 . An important development in this context was made by Kanzieper u&, where the link between 
nonlinear replica a-models and Painleve hierarchy of exactly solvable Toda lattice equations^ was established. More 
precisely it has been shown, that the replica partition function of the Gaussian unitary ensemble (GUE), which in it's 
original and dual (in the replica space) cr-model^ representation satisfies the Toda lattice equation (TLE), can be 
reduced to Painleve transcendent tpi\£. From here exact nonperturbative results were obtained demonstrating the 
ability of the method to overcome the critique of replica trick 8 . In the chiral unitary ensemble (chUE) , the replica 
trick on the basis of TLE was analyzed inML. 

In this paper we show how the above picture may arise and work in a unitary ensemble defined by probability 
distribution function P(H) = exp[— trV(H)] with the confining potential 



that behaves as V(x) oc ln 2 (|x|) at |x| — > oo. Here x — sinh/?x/2 and the parameter q — e - ^, with (3 > 0. 

Defined and solved using the orthogonal polynomials method inifi, this unitary random matrix model has attracted 
considerable attention during the last decade. One reason for that is that it leads to the eigenvalue distribution which 
for not very large (3 < ir 2 is intermediate between the Wigner-Dyson — and the Poisson distribution and has the 
features typical of the critical level statistics ^Li£. 

Another reason, perhaps more important, is that the potential Eq.lQJ fits very well the function a[ln(l + bx 2 )] 2 in 
a wide region of the variable x. This behavior was obtained numerically in the tight binding Anderson model with 
random site energies as a distribution potential of variables Xi, i = l,...,N, connected with the conductance as 
g = + x 2 ) _1 — , reflecting a log-normal behavior of the conductance g in a product of random matrices—. 

Therefore, by treating matrices H as H 2 = (TT + + T + T — 21) /4, where T is the transfer matrix characterizing the 
disordered conductor, the model with the potential can be considered as a qualitatively correct phenomenological 
model for the Anderson insulator with the parameter playing a role of the ratio of the system size to the localization 
length L/£ and x/£ playing a role of the Lyapunov exponent (J is an identity matrix). At large (3 the probability 
density p(x) is strongly peaked near positive integers exhibiting the phenomenon which is known in theory of 
quasi-ld Anderson localization as transmission eigenvalue crystallizationi^^^^. 

Up to now the model Eq.Q was investigated only within the approach of orthogonal polynomials which was 
successful because the corresponding set of orthogonal polynomials (q-deformed Hermite polynomials) is explicitly 
knowni&. Though this method gives an exact solution to the problem it does not tell us anything about the connection 
of this problem to exactly integrable lattice models. On the other hand, such a connection is expected by analogy 
with the Wigner-Dyson RMT which is solvable by the ordinary Hermite polynomials. Thus the main question we 
ask is as to whether or not the RMT Eq. is related with exactly integrable lattice models and what precisely these 
models are. 

As a first step we compute the replica partition function exactly using the pole structure of the distribution function 
P{H). Indeed, the form of the potential allows us to integrate over poles in complex plane and obtain the new exact 
expression for the partition function as an infinite sum of residues. Such a result of nonperturbative nature turns out 
to be exactly deducible from the dual representation of the partition function in a form of the Itzikson-Zubcr (IZ) type 

matrix integral over homogeneous factor space of pseudounitary supergroups g^^M^JVM) (Stiefel manifold) with 
M — > oo . Thus one of the results of the paper is a calculation of this specific IZ integral by use of Duistermaat-Heckman 
theorem. To the best of our knowledge this integral was not calculated earlier in the literature. 

One of the main results of the present paper is the establishment of the link of RMT Eq. (JJJ with integrable lattice 
models. For the original replica partition functions we infer the same Toda lattice equations as in the case of the 
Gaussian RMT 2 . However the dual partition functions Z n ^ obtained by a super-symmetric mapping onto a variant 
of the replica a model (color-flavor transformation^) has a more delicate relationship with the Toda lattice models. 
Namely, semi-infinite hierarchies of fermionic and bosonic flavors involved in this mapping manifest themselves as an 
extension of graded TLE first found in 21 in the context of chiral unitary ensemble relevant for QCD. 

Finally, we analyze the new representation for the two-level correlation function R2{x) and the probability density 
p(x) an d discuss the phenomenon of ^-crystallization. 

The paper is organized as follows. In the second section we make a color-flavor transformation in the model 
under consideration and give its sigma-model formulation defined on the homogeneous factor space of pseudo-unitary 
supergroups. We calculate the replica partition function in both, color and flavor spaces respectively. In the third 
section we infer Toda lattice equations in both spaces. In the fourth section, by use of the method of pole integration, 
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(1) 



n=0 



we calculate exactly the density of states and the two point correlation function at large N. Comparison with known 
results follows. Section 5 contains our conclusions. The calculation of the Itzikson-Zuber matrix integral over Stiefel 
manifold and the derivation of TLE in flavor space are presented in two Appendices. 



II. THE REPLICA PARTITION FUNCTION. 



The replica partition function of the model Eq.JU, defined by Z n ^(e) = (det™(e — H))h, can be represented as an 

1 T~\f CKTT"C1 1 



integral 

M 

{C M ) N / dHdet n (e - H ) TT — -j- f — — T , (2) 

where fi k = cosh^, C M = (Il£Li V) -1 and the limit M 

— > oo is understood. The main idea of the replica trick is 
the possibility of analytic continuation of Z n ^ to non- integer values of n, which was shown to be exact in GUE 2 . 

Diagonalizing Hermitian NxN matrices H = V diag{Ei] and taking into account the Vandermonde determinant 
one can calculate Eq.JSJ by closing the contour of integration at infinity (which is correct for n < 2M — N) and using 
the Cauchy formula of pole integration: 

Z n>N (e) = {2tt) n (C m ) N ^n^^m-n!^ Z> ( 3 ) 

where the sum is taken over all possible subsets {k\, . . . , k^} of the set of indices k 6 {1, . . . , M} of the poles fi k + iEi 
on the upper half plane and {pi, . . . ,P2M-n} = {1, • • • , M, M + 1, . . . 2M} \ {fci, . . . , kjy} includes the numeration of 
poles /ife — iEi with the condition Hm+o — ^Hj, 3 = 1> • ■ -M. This formula will be a starting point of the further 
analyzes. 



A. Color—flavor transformation 



Now we define another, dual representation of the partition function in which the replica space and the space of 
matrix indices interchange. Usually (see e.g^) such a representation is obtained from the fermionic replica sigma- 
model derived using the Hubbard-Stratonovich transformation. 

In our case the potential is not Gaussian and the Hubbard-Stratonovich transformation is not useful. But 
nevertheless it appears that one can make a transformation of the partition function given by an integral over original 
NxN Hermitian matrices (following^ let us call the space where they are acting a color space) and pass to the 
integration over matrices, which are acting in replica space (let us call it flavor space). Similar type of transformations 
effective in supersymmetric er-models^ are called color-flavor transformations. 

In order to make the color-flavor transformation let us transform the determinants in the integrand of Eq. J5J into 
the Gaussian integral over superfields $° = ({*"}; {^i^ ,k }; {4>T" >,P }) with a = {a; k;p} = {1 . . . n; 1 . . . M; 
1 . . . M}, where n is the number of fermionic components while M is the number of bosonic ones 4>f^ ,k , ^ '. The 
supersymmetric structure of the fields <I>° which can be regarded as i — 1, . . . N vectors in (n+2M) dimensional flavor 
superspace, is dictated by the fact that the partition function Eq. Q is a product of a certain number of determinants 
and inverse determinants. Using these superfields we obtain 

= \ 27r) 2MN / dHRd*t a d*S<*p(-*t a Vv*i)' (4) 

^ * i,a 

The indices i, j — 1 ... N correspond to the original color space while indices a and b are forming a basis in the flavor 
space. The matrix V°j b = diag{[ieSij — iHij]d a p; [Hk&ij + iHij]5k p ; [^k^ij — iHij]5k p } is diagonal with a,j3 = I . . . n; 
k.p = 1 . . . M and Kronecker <5-s. As one can see the integral over Hermitian random matrices H in the expression 
® simply produces a Dirac delta-function {2tt) n ^ n+1 ^ 2 H id 5[<S>+ a g ab <i> b j } with g ab = diag{S lj: S kp: -S kp }. 

For further treatment it is convenient to make use of the following trick: we introduce a set of real numbers 
wi,...,(Jjf and instead of the delta-function FJ i . 5[§f a g ab § b ] we substitute lim^.^o Hi j S[^f a g ab ^> b - LUiUi^Sij}. 



Then after the rescaling of the superfields $f — ► one obtains 



■ J] mt a 9 ab <S>) ~ Sij] eM-^M a 9 ab V^ C j), (5) 

with 

/ ieSijS a \ 

(6) 

and the summation is understood over repeating indices b. The argument of the new S- function in Eq.0 imposes 
that all possible field configurations contributing to the integral are N normalized vectors in the (n + 2M) dimensional 
complex superspace with metric g ab . By definition this space can be called complex Stiefel manifold and is equivalent 
to the homogeneous factor space S = su^u\m-n\i) °^ pseudounitary supergroups. Therefore the formula (J^J is 
nothing but a replica a model formulation of the model 10 in the form of graded IZ type integral over Stiefel manifold 
in the (n + 2M) dimensional flavor superspace: 




- (2 JmI C ^ )/2 afi^ I dUeMStr(aU-^VU)}, (7) 
where U^ 1 = gU + g, the matrices V — diag{iei . . . ie n , fii . . . [Im, — Ml • ■ ■ — Mm}; = diag{0 . . .0, w\ . . . ujJj, ... 0} 



2M—N 



and the limit ei, . . . , e„ — > e is supposed. 

The evaluation of standard (super) IZ integral over the group manifold is based on the analyze of heat kernel equation 
in the Hermitian matrix space and on calculation of induced metric on the subspace of diagonal matrices^^jSt. j n 
our case, since the tangent to the Stiefel manifold vector defined by dH = UdU^ 1 \ GdG^ 1 , with U G SU(n \M,M), 
G G SU (n | M — N, M), the corresponding induced metric does not coincide with the Vandermonde determinant, but 
with some modification. 

Another elegant calculation of IZ integral over unitary supergroups has been made in2£ by use of Duistermaat- 
Heckman theorem. In^ it has been shown, that the theorem can be applied to pseudounitary supergroups as well. 
In Appendix A we apply this theorem to our IZ integral over manifold SU(n | M, M)/SU(n | M — N, M). Leaving 
the mathematical details of calculations aside we present here only the answer. Let us introduce a family of maps 
a : {1, 2, . . . , N} — * {1, . . . , M} of the indices of the original color space into the subset of indices of the positive 
root part of the flavor superspace. Correspondingly a will denote the complement part of the a: a = {1, . . . , 2M} \ 
{cr(l), . . . , cr(N)}. Then the answer is (see Appendix A for details) 

/ dUeMStrinU-'VU)} - ( 27r )2^-^(iV-i)/2 y d& f Xp ~f k ^ ] (8) 
Jues f£ A CT (ft)A CT (lO 



where a, ft = l,n; k,p — 1, N are fermionic and bosonic components correspondingly, the sum is taken over all possible 
maps a and A ff (w), A CT (T^) are defined by the induced metric as 

T-riV T-f2M-N, w \ 

a t \ U»=illg=i v^W ~ v s{p)){v<T(i) - , . 

A<j(^) = T7 . (9) 

n^iiiwK-^i)) 

The formulae ||HJ| and @ present the extension of IZ matrix integral for the case of integration over Stiefel manifold 
S. It is clear, that the expressions (|SJ 10 were obtained for integer values of n. Now taking the limit uj i — > in (J7J) 
we will reproduce the same formula 10 for replica partition function as it was obtained in original color formulation. 
This demonstrates the equivalence of the color and flavor formulations of the model under consideration. 



III. TODA LATTICE EQUATIONS 



In the present section we derive Toda lattice equations in both, color and flavor spaces respectively. The first part is 
devoted to the so called r function, defined for a unitary ensemble with any given distribution function, and TLE for it 
in color space with initial conditions corresponding to the RMT Eq. 0. In the second part we show how the fermionic 
and bosonic flavors involved in the replica partition functions Z n ^, connect them with integrable hierarchies. 



A. TLE in color space 



It is easy to see, that the replica partition function J5J is the T — > limit of so called r-function t^(T, e) — 
J dH expTr[TH + V(H)j of the model with the potential V(H) = nln(e — H) + V(H), which can be regarded also 
as a matrix Fourier transform. N x N matrix T, being conjugate variable to energy matrix H 7 can be regarded as an 
effective time. The r-function can be calculated exactly as2£ 



tn{T; e) 



A(T) 



det[Ws(tj)]5,j=o,...,N-i, 



where A(T) is the Vandermonde determinant corresponding to the matrix T and 

W s (t)= [ d\\ s exp{\t + V(X)} 



(10) 



(11) 



with 5 = 0, . . . , N - 1. In the U -> limit (JTOJl gives us the replicated partition function Z n! jv(e) = Tiv(0, e) and 
T/v(i, e) = tn(T, e) \{ti}=t can be expressed via Hankel determinant 



N ~ 1 xs 

r N (t,e)= det [d> +s W (t)]. 

s,j=o 



Then it is clear— (Si, that Tjv(t, e) will satisfy the TLE 



T N (t,e)d^T N (t,e) - (d t T N (t,e)) 2 = T N+ i(t,e)T N -i(t,e) 



with the initial condition 



(12) 



2M 2M 

n (0,e) = 27rC(M)^[(e-^ fc )7 J] (j*k ~ lh) 

k—1 p—l;p^k 



(13) 



B. TLE in flavor (supersymmetric) space 



After some algebraic manipulations the replicated partition function @ can be rewritten as 

(C M ) N .. 1 



,N 



(e) 



lim 



V 2 m(^i ■ ■ ■ M2Af) e =^ e V„(ei . . . e n ) 

/A~l) ^JV+n(ei ■ ■ - e n,/V(l) ■ ■ ■ Va(N))V2M~N(Llg(l) ■ ■ ■ Vs(2M-N) ) 



M 



(14) 



where V({. . . }) is the Vandermonde determinant defined by the set of arguments {. . . } and ir(cr) is the parity of 
permutation of 1, . . . , N, N+ 1, . . . , M to cr(l), . . . , a(N) 7 a(l) . . . , a(M — N). This form of Z n ^{n) can be understood 
as a formula for a determinant calculated by minors 



(*-™C M ) 



JV 



lim 



I- 1 det [A] 



V2m(mi • ■ ■ M2m) V„(ei . . . e„) 
where we have introduced the (n + 2M) x (n + 2M)matrix A as 



(15) 



A = 



( h Sih 

h S 2 h 



In 3n In 

h 6J1 



Im SmIm 




rn+N-1 



h 



u n 1 n 








h 



() 1 





2A/-7V-1 



h 



V o 



Im+i 



12M 



J2M-N-X T 

M Im 

2M-N-1 f 
M+l J M+1 



F2M-JV— 1 
J 2M 



(16) 



12M 



involving the functions I a — exp(ie a ), a = 1, ...n; I r = exp/i r , r = 1, . . . 2M; / = HaV=i ^alr and the notations 
*a = £a^; 5 r = Mr Define <5 1+2M = e£ + X^fr). 

Now, following the technique developed in2i, the modified graded Toda lattice equation for replicated partition 
function Z n ^(e) will follow directly from its expression (|15|l as a determinant and two sets of Sylvester identity 2 ^ for 
the matrix A (see Appendix B for a details). As a result, the following equation is established 

n-id t 8 1+2M HZ n . N {e)] = 

where n < 2M — N, when Z n ^ is finite. It has a graded form since contains derivatives over both, fermionic and 
bosonic parameters. One can see here, that in the replica limit n — > the right hand side (RHS) of this equation 
involves a factorized product of fermionic and bosonic partition functions Z n= \ and Z„ = _i correspondingly. Moreover, 
the second term in RHS shows that indeed fermionic, bosonic and supersymmetric partition functions for any N belongs 
to the same integrable hierarchy. 

In the limit N — > oo (but N < M), when Z„jv-i = Z„jv+i = Z riiOC , the equation ljl7Jl simplifies to 

■a rl + OO 1 Try I \1 r, Z n + 1 .OO ( e ) Z n — \ oo (g) f-\o\ 

n ~ id t S ]n[Z nt00 (e)\ = 2n — — — , (18) 

where 5 1+ ~ = + E"i(^r). 

Similar to 1)18(1 type of equations, but for a finite amount of mass parameters [i were considered earlier inSH^jSi^. 
The factorization properties of Toda lattice equations in the replica limit have been analyzed inSi. 

As we see the equation (|17|) differs essentially from the known equation JTSJ, which defines integrable Toda lattice 
modeli Since the equation (|1 7|> appears for the exactly calculable partition function J2J of the unitary matrix model 
with potential lfT)l. one can conjecture the existence of corresponding integrable lattice chain model. It is definitely 
interesting to investigate the interplay between color and flavor degrees of freedom in the equation l|17fl further and 
analyze its consequences. 



IV. DENSITY OF STATES, MEAN CONDUCTANCE AND TWO POINT CORRELATION FUNCTION 

In the original articles^^ where the model with the potential JQ) was initiated, the authors, by use of q-Hermite 
orthogonal polynomials, have analytically calculated the N — > oo asymptotic of the two point correlation function. 
For the appropriately renormalized two point kernel in the bulk of the spectrum exact expression approximately 
becomes (except very close to the origin) 

"Ml- 27rsinh[(£-7,)/3/2]' ^ 

which is translational invariant. In its /3 — *• limit K(^,r]) reproduces the universal behavior of Wigncr-Dyson RMT, 
while for the large (3, as authors ofiS have argued, the model approaches to the uncorrelated systems with Poisson 
distribution. 

In our approach the calculation of the density of states p{x) will be reduced to Cauchy integration over poles of the 
distribution function in the 

i r°° Ntr ) N r N N M ' N 1 

ZoM e ) J-oo Z . N (e) J A± 1^ i^^^ k + lx .^ k _ lx .} 

where the normalization factor Zo,jv(e) is defined by the expression Q, for the case n = 0. In order to analyze the 
large N limit we will consider N — M — > oo. In this case the calculation of the Cauchy integrals over poles in our 
problem essentially simplifies. From © one can deduce 

M 

ZoMe) = {^) M {C M ) M M\ TT -. (21) 

kP ii iMKi+ihi) 



Now, in order to take the Cauchy integral over poles in l(2Q|l we close the integration contours at infinity of upper 
half plane and notice, that because of presence of Vandermonde determinant, the integration variabes Xi, i = 2, . . . , M 
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FIG. 1: x eigenvalue densities at /3 — 6 (dashed line) and j3 — 20 (solid line). The connection with transmission eigenvalues of 
T is given by T = cosh[/3x/2]~ 2 



will utilize M — 1 poles, all different. Only one pole will be left untouched and we will have the summation over it in 
the result: 

Constituents of this formula have the following meaning. The factor (x 2 +fi1) Yl q ^(x+ifJ, q ) in the denominator comes 
from the distribution function of the remnant eigenvalue x, the part of which, namely Oq^fc — was cancelled 
with the similar bracket in the Vandermonde determinant. The factor Hg^fc( :E — m the nominator is the remnant 
from the Vandermonde determinant. The factor TT^Si. ^ k+tlq \ is left after mutual cancellations between the residues 
of integrated poles, the Vandermonde determinant and the normalization factor Zo,./v(e) defined by formula (|21(l . The 
second part of the formula (122[l is an identity. 
Hence, we have obtained a simple formula: 

1 M 1 

^) = ^X>(^)- (23) 

One will arrive at the same formula (|23J) taking a replica limit for the one point Green's function G(x) — 
lhn n ^on~ 1 dZ Tlt M(x)/dx and calculating the DOS. 

The model under consideration with potential Q first was considered inifi, in order to describe the disordered 
conductors connecting matrices H with the ensemble of random transfer matrices T as 

H 2 = i(TT+ +T+T - 21). (24) 

Therefore substituting x — sinh[/3x/2] (eP x ^ 2 is the eigenvalue of T) into the p(x) one will obtain an exact expression 
for the density of x' 

P(X) = ^P(x). (25) 
ax 

which exhibits crystallization at the integer values n, as it is presented on Fig.l for two values of f3. This phenomenon 
for different ensembles was observed earlier in the articlesi^r^^i. 

The eigenvalues X4 of H and Xi are connected with the conductance g = tr(tt + )/2 (t is the transmission matrix) as 



JV N 

i . , ^ cosh[/3 Xl /2] 
and by the same type of Cauchy integration over poles as in p(x) one can find the average conductance 



9 ^1 + x 2 £ cosh^W^ 2 (26) 



M M ^ 

^ = E TT7Z = ^l + coshfn/3/21 2 ' ^ 

n=l r " n=l 



R 2 [f, r?] 



l 




FIG. 2: The two point correlation function i?2 ~ ty] calculated exactly by Cauchy integration over poles and after unfolding 
with p(x) (dots). The solid line represents i?2 — rj] obtained from the formula I19L both at (3 = 0.3. 



In the same way one can calculate the two point correlation function R2[x, y], again in the limit N — M — > oo 

R 2 (x,y) = —± 7 - [ dHtr6(x - H)trS(y - H)P(H). (28) 

By Cauchy integration over poles, in this case, we should now leave untouched two poles of the distribution function 
P{H) (see formula (J2UJ). Therefore, in the result we will have a sum over all possible positions m and n of that poles. 
The answer reads: 



R 2 (x,y) = ^J2 



M 

Prn Pn 



n 



^ ^ + Vl)( x2 + tfn){y 2 + ^IW + Mm) 

ipk)(p m +fJ>k)(y - «Mfc)(Mn + Mfc) 



(29) 



Details of the calculation of R2(x,y) repeat the detales of the one for the density of states p(x), presented above. 
Since here we hold two eigenvalues x and y untouched, the formula l|29() is similar (but not equal) to the product of 
two structures 112211 of p(x) with corresponding adjustment of the summation procedure. Unlike to the case of p{x), I 
was not able to find a simplification of the final formula (|29(l , as it was done in the second part of Eq. (1221 . 

In spite of explicit presence of the imaginary unite i in the expression l)29|l the resultant R2(x,y) is real. After 
unfolding, namely after the re-scaling x — » £, defined by d£/dx — p(x), i?2(£> v) is pictured in Fig. 2 by dots for (3 = 0.3. 
Moreover, the numerical analyze of the expression Ij29(l after unfolding shows that i?2 (£,??) in fact is a function of 
difference £ — rj. 

Numerical analyses show, that consideration of other potentials, say by taking the product in formula J2J from 
k = or k = 2, after appropriate unfolding gives us the same behavior as approximately defined by formula (|19fl . 



Conclusions 



We give a sigma- model formulation of the RMT with confining potential Q . The replica partition function of the 
model has been calculated exactly in both, original color space and in flavor space (after a color-flavor transformation) 
with the same result. In both spaces we have established connections between the RMT Eq.Q and exactly solvable 
lattice equations. A generalization of graded TLE is defined and solved. 

We have exactly calculated the probability density and the two point correlation function of the model at N — > oo 
and confirmed the known results obtained earlier. The effect of crystallization of transmission eigenvalues has been 
observed. 
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Appendix A. Evaluation of Supersymmetric Itzikson-Zuber Integral over SU(n \ M, M) / SU(n | M — N,M) 



In this appendix we present a simple derivation of the supersymmetric extension of the Itzikson-Zuber formula over 
Stiefel manifold SU(n \ M, M)/SU(n \ M — N, M). 

l(X,Y;r)= [ dU exp{r StriXU^YU)} (30) 

JUeSU(n\M,M)/SU(n\M-N,M) 

In this calculation one may follow the standard technique 2 ^ developed m 2 ^^ for unitary supergroups based on the 
solution of the heat kernel equation. In^ 5 . it was argued, that the heat kernel equation technique can be applied to 
pseudounitary supergroups. However we will present here the calculations based on the application of Duistermaat- 
Heckman theorem to integrals over supermanifolds, which states that under certain conditions the quasi-classical 
saddle-point approximation is exact^. This theorem was successfully applied to IZ integral over unitary supergroups 
U(n | m) inS and pseudounitary supergroups U(n \ m, k) irM. In what follows we will assume, that this criteria are 
met by our integral (|38|l as well and evaluate the integral in the saddle-point approximation. The fact, that final result 
reproduces the same answer (@J for Z n ^(e)^ as we got by simple calculation of the integral J3J by the Cauchy formula 
confirms, that this assumption is correct. The same answer one will obtain following the standard techniqu e] 25 ) 26 . 

We need to find the extremum of the function 

A(U) = Str^XU^YU), (31) 

where U € SU (n \ M, M) / SU (n | M — N, M) and X and Y are diagonal matrices from the expression JJJ). It is easy 
to find the saddle-point equation 

[X, U^YU] = (32) 

from the condition SA{U) = Str{{XU~ l YU - U- 1 YUX)U~ 1 SU} = 0. 
Since we have X — diagfO^O, x\ . . . xn,Q_^0} and 

n 2M—N 

U € SU(n | M,M)/SU(n \ M - N, M) the only solutions U of © (up to irrelevant U € SU(n \ M - N, M)) are 
those matrices, which permute the diagonal elements of the matrix Y in three- fermionic, positive and negative metric 
parts of the bosonic sectors separately. In other words the set of solutions of (|3*2*|l is 




, n 



U^=\ U kp I, k,p=l,...,M (33) 



ILfi \ a,/3 = l, 

| II/;. 

' ' k', P '=l,...,M 

where Jl Q/3 and Ilfcpdlfe'p') are the elements of the permutation groups S n and Sm of the fermionic and bosonic 
sectors respectively. Notice, that the set of permutation elements i, j — 1, . . . N coincides with the set of maps 
a : {1, . . . N} -> {1, . . . M} defined earlier. 

We now set U = U^e tL in Ij31|l . with an infinitesimal Hermitian supermatrix L = U~ 1 5U G su(n \ M, M) \ su(n \ 
M — N, M) and expand the function A{U) up to the quadratic order in L. In this factor space L has only following 
nonzero matrix elements Ly, i, j — 1, ... N; Lik, i — 1, . . . N; h = 1, . . . 2M — N and Li ai i = 1, . . . N; a = 1, . . . n. 
In the second order over L, A(U) is 

-4(17) = StrXU {0) YU^ + StrXLU {0) YU {0) L- Str^X.U^YU^L 2 

= StrXU {0) YU {0) + ^Str[X,L][U {0) YU {0 \L]. (34) 

Due to Duistermaat-Heckman theorem in evaluation of IZ integral (l.'U)ll we should take the sum over all extrema, 
therefore substituting Ij34(l to (|30J) we obtain 



N 



J(X,Y;r) = ^2 ex P( r ^2 x iy<r(t)) / Yl dLi 3 II dLik II dL 



N 



N 



N.2M-N 



,k=l 



i,a — l 



x exp r 



N.2M-N 
i,fc=l 



(35) 



where y Q are fermionic components of the diagonal matrix Y and a is the complement to a permutation, defined after 
the formula Q. Here we essentially have used diag{(L^J), X\ . . . xat, (T^T)} form of the matrix X . 

n 2M-N 

Evaluating now the Gaussian integrals in (|35|) over complex bosonic Ly, i ^ j\ Li_k, i k and the complex 
Grassmann variables Li a , we arrive at 



T(X,Y;r) = (2 7T ) 2MN - N( - N - i y 2 (~r) Nn ~ 2MN+( - N ~ 1 ' )N/2 
det[exp rxiy^)] 



A a (X)A a (Y) 



(36) 



with 



T-riV TT 2M-iV/ \ / \ 

rui n P =i (Mj) - ^(p))(^ g (i) - ^(j)) 
n^UiiiiliK-'Mi)) 



i/ = x, y. 



(37) 



Appendix B. Graded Toda Lattice Equation 



In this Appendix we show how the graded Toda lattice equation for the replica partition function Z n n(c) appears. 
First we would like to demonstrate, that Z n ^(e) is a ratio of certain determinants. The expression can be regarded 
as e a — > e; a = 1, ... to limit of n different fermionic eigenvalues. For the convenience and compactness of formulas we 



will drop from further transformations the constant factor (2n) (Cm) 



\N-nN 



and introduce the notation e a = 



for a while. Also, instead of fixed [i\ . . . [iMi — A*i, ■ • • — [Im set of parameters of the model, we consider here general 
situation with different ji}., k = 1, . . . 2 M. Condition fik+M = — Mfe can be easily imposed at the end. We have 



nLiIl^iK - /Mi)) 



'"•" <e) = w-^n£ I nS-"( 



{ff} lli=lllj=l (Mo-(i) ^s-(j) ) 



Um „„ . . —jv -rr2A/- 

{e ° } ^ e w ria</3=i(e Q - e/j) n<=i n J= i 



(38) 



where the nominator and denominator have been multiplied by the same expressions. Introduce now a bosonic 
Vandermonde determinant: 



V s ({xx,...x s }) = Y[ (xi - Xj) = Det 

i<j=i 



1 X\ 

1 x 2 
1 I, 



(39) 



The expression (|38|) for the partition function can be rewritten as 

. V„ + at({#i, . . . n+N })V2M-N({^a(l), ■ ■ ■ Ha(2M-N)}) 



Z n<N (e)= lim ]T(-1) J 



V„({ei, . . .e„})V 2 Af({Mi, ■ ■ -M2m}) 



(40) 



where {#1, . . . 9 n +N} — {ex, . . . e n , . . . /^(jv)} an d ( — 1) P i s the parity of permutation 

P : {e\, . . . e n , Hi, . . . ^2m} — ► {ei, . . . e n , . . . fi a (N), Mff(i): ■ ■ ■ Ms-(2M-A')}- (41) 

The denominator in the expression 1|40[1 does not depend on a and can be taken out of summation, while of 
nominators is nothing but the determinant of the matrix Aq : 



A 



/I ei 

1 e 2 

1 e„ 

1 m 



e^- 1 
e* +N - x 



Mi 



mi 



1 o 
1 i 





Ml 



1 Mm Mm 




V 



Mm^ 1 1 







Hm 

1 A^M+1 
1 M2Af 



1M-N-1 
Ml 



2M-W-1 
MAf 

2M-JV-1 
Mm+1 



2M-JV-1 / 
M 2 Af 7 



(42) 



calculated by N + n and 2AI — N dimensional minors and by use of the formula (|39J) . 
Now let us consider the operator 5 X — x-^ introduced earlier and observe that 

e- x (5 x ) p e x = x p + V {p - 1] {x) (43) 

where V^ -1 '^) is the polynomial of p — 1 order of x. From this it follows, that if we replace the monomials x p in 
the matrix Aq with e~ x {5 x ) p e x , we will obtain the matrix A defined in Their determinants relate as 

det[A ] = det^A] (44) 

and define the nominators in the expressions Ijl5(l and H40(l for the partition function Z n ^(t). The notation / = 

ex P(ELi e s + J2k=i Mfe) is defined just after formula l|15|l . 

The next step is the consideration of the degenerate limit e a — > e; a — 1, ... ,n. We put e\ = e and consider Taylor 

expansions f(e a ) — /(e) + (e Q — e)<9 e /(e) + . . . ^ a (n-\)\ — ^" _1 /( e ) + • ■ ■ around e of all entities in the upper-left n x n 
corner of the matrix A in Ijl5(l . It is easy to see that in the limit e a — > e; a = 1, . . . ,n, higher than n. — 1 orders 
of the expansions will not contribute into the matrix A, while the coefficients (e a — e) k will form a Vandermondc 
determinant (|39[) . which will cancel the same V n ({ei, . . . , e„}) in the denominator of Z n ^(e) in the formula (|4UI) . 
Therefore we obtain 



Zn,N{t) 

with / = exp(ne + ^ r /^ r ) and B n ^ equal to 



V 2 m({mi, • ■ -M2m}) 



( 



h 



5 n+N-lj i 



h Sji 



(n-l)\ U e °1 



-lxn+N-1 



n+N-1 



h 



Im 




$m!m 




cn+N-l f 








/i o 

h 

Jm 
Im+i 

I2M 



(45) 



72M-N-1 f 



Af+1 J A^+1 



?2M-JV-1 ; 
°2M 



(46) 



'A/ 



In this expression I\ — exp(e) and 8\ = ed e . 

In order to derive the Toda lattice equation (|17|l we will follow the technique presented in the article^! in connection 
with QCD partition function. The method is based on use of the Sylvester identity^. For the determinant of any 
given matrix D 

Ci,jCp,q Ci,qCp.j — det[D]CiJ-p,q, (47) 

where Cjj is the cofactor of the matrix element ij: 

ddet[D] 

a <> ~ ~dD — > (48) 
and C'ij-p^ q is the double cofactor of the matrix elements ij and pq: 

d 2 det[D] 

™ <>!),. ,i)l) p .,- l4Jj 

Let us consider now the (2M + n + l) x (2M + n + l) dimensional matrices i? n +i.jv and B n+ i t N-i and write Sylvester 
identities (|47|l for different set of indices i,j;p,q for each of them. 

For the B n+ i t N we take i = n + 1, j = n + N + 1, p — n, q = n + N and obtain 

.n+N 

= det{B n+ i t N]C n +l,n+N+l;n,n+N ■ (50) 

From analogous to (|46|l structure of B n+ i^ follows that 

Cn+l.n+N+l = det[B n at], Cn.n+N+1 = —d e det[B n *f], 

n 

C n +l,n+N+l;n,n+N = det [B n -i t ff] , (51) 

C n ,n+N = —9 e det[B' N ], C n +i ,n+N = det[B' n N ], 

n 

where the matrix B' n N is defined by acting by the operator 5 1+2M — e + X^r=i ^ r sf~ 011 ^ ne ( n + column of 

the matrix B n> ^. 

Now consider the Sylvester identity for the matrix B n+ i^-i- Setting i = n+1, j = n + 2M +1, p = n, q = n + 2M, 
we will have 

/^if 

u n+l,n+2M+l L/ n,n+2M — Ly n,n+2A/+l°ri+l : nH-2M 

= det[B n+ i^ N -i\C' n+ln+2 M+l;n,n+2M- 

(52) 

As in the previous case and from the structure of j^l for B n+ \^-i it follows that 

Cn+l,n+2M+l = det[B n ^], ^'n,n+2M+l = ~&edet [B n> If] , 

C n +l,n+2M+l;n,n+N = [B n -l,N+l] , (53) 
C n ,n+N ~ —dedet{B" N ], C' n+l ^ n+N = det[B' r ' l N ], 

where the matrix B'^ N is defined by acting by the operator S 1+2M on the (n + 2M)-th column of the matrix B n ^. 
It is easy to observe, that 

det[B'^ N ] + det[B^ N ] = 5 1+2M det[B n , N }. (54) 

Inserting now expressons H51(l and l)53|l into the Sylvester identities (|5()[) and Ij52() respectively taking their sum and 
using |[5Hl we get 

det[B ntN ] 2 ^5 1+2M \n[det[B ntN ]] = (55) 
oe 

= n det[B n+liN ]det[B n - ltN ] + n det[B n+ i tN -i]det[B n _x !N+1 \. 



Finally, using the relation l|45|l between the partition function and the matrix B n ^ we obtain the desired graded 
Toda lattice equation as 

n + d e 6 1+M \n[Z n , N {e)] = 

^n+l,Af(e)^n-l,Af(e) N-l (g)^n- 1, N+l (g) 

which coincides with l|17(l after substitution e = ie. Notice that the Vandermonde determinants, V^mCImii ■ • -M2m}), 
do not contribute to the left hand side of the Toda lattice equation since S 1+2M ln[V2M • • ■ M2m})] is a number, 
while the contribution of I = exp(rte + X)r=i A 4 * - ) gives addition of n. 

If we take N — M, then, since Bn-^jv+i = 0, the graded Toda equation (|56|l simplifies. 



1 M. L. Mehta, Random Matrices (Academic Press, Boston, 1991). 

2 E. Kanzieper, Phys. Rev. Lett. 89, 250201 (2002). 

3 M. Toda, J. Phys. Soc. Jpn. 22, 431 (1967). 

4 S. F. Edwards, P. W. Anderson, J. Phys. F 5, 965 (1975). 

5 F. J. Wegner, Z. Phys. B 35, 207 (1979); 

6 K. B. Efetov, A. I. Larkin, and D. E. Khmelnitskii, Zh. Eksp. Teor. Fiz. 79, 1120 (1980). 

7 P.J. Forrester and N.S. Witte, Commun. Math. Phys. 219, 357 (2001). 

J. J. M. Verbaarschot and M. R. Zirnbauer, J. Phys. A. 18, 1093 (1985); M. R. Zirnbauer, cond-math/9903338 
K. Splittorff and J. J. M. Verbaarschot, Phys. Rev. Lett. 90, 041601 (2003). 

K. A. Muttalib, Y. Chen, M. E. H. Ismail and V. N. Nicopoulos, Phys. Rev. Lett. 71, 471 (1993). 
Y. Chen, M. E. H. Ismail and K. A. Muttalib, J. Phys Condens Matter 4, L417 (1992). 
Y. Imry, Europhys. Lett. 1, 249 (1986). 
K. A. Muttalib, Phys. Rev. Lett. 65, 745 (1990). 

J. L. Pichard, N. Zanon, Y. Imry, and A. D. Stone, J. Phys. (France) 51, 587 (1990). 
E. Bogomolny, O. Bohigas and M. P. Pato, Phys. Rev. E 55, 6707 (1997). 

16 M. E. H. Ismail and D. R. Masson, Trans. Am. Math. Soc. 346, 63 (1994). 

17 V. E. Kravtsov and K. A. Muttalib, Phys. Rev. Lett. 79, 1913 (1997). 

18 V. E. Kravtsov and A. M. Tsvelik, Phys. Rev. B. 62, 9888 (2000). 

19 M. R. Zirnbauer, J. Phys. A. 29, 7113 (1996). 

20 K. Splittorff and J. J. M. Verbaarschot, Nucl. Phys. B 683 467 (2004). 

21 K. Splittorff and J. J. M. Verbaarschot, Nucl. Phys. B 695 84 (2004). 

22 C. W. J. Beenakker and B. Rejaei, Phys. Rev. Lett. 71, 3689 (1993); 
K. Frahm, Phys. Rev. Lett. 74, 4706 (1995); 

C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 

23 A. Kamenev and M. Mezard, J. Phys. A 32, 4373 (1999). 

24 C. Itzykson, J. B. Zuber, J. Math. Phys. 21, 411 (1980). 

25 T. Guhr, J. Math. Phys. 32, 336 (1991). 

" J. Alfaro, R. Medina and L. F. Urrutia, J. Math. Phys. 36, 3085 (1995). 

A. Morozov, Usp. Fiz. Nauk 164, 3 (1994) [Phys. Usp. 37, 1 (1994)]; R. Dijkgraaf, [arXiv:hep-tn /9201003 . 
J. J. Sylvester, Compt. Rend. Acad. Sc. 54, 129 (1862). 

G. Akemann, D. Dalmazi, P. H. Damgaard and J. J. M. Verbaarschot, Nucl. Phys. B 601, 77 (2001). 
S. Kharchev A. Marshakov, A. Mironov, A. Morozov, A. Zabrodin, Nucl. Phys. B 380, 181 (1992). 
A. Lamacraft, B. D. Simons, M. R. Zirnbauer, Phys. Rev. B 70, 075412 (2004). 
A. Schwarz and O. V. Zaboronsky, Commun. Math. Phys. 183,463 (1997). 
R. J. Szabo, Nucl.Phys. B 598, 309 (2001). 
Y. V. Fyodorov, E. Strahov, Nucl.Phys. B 630, 453 (2002). 
A. V. Andreev, B. D. Simons, Phys. Rev. Lett. 75, 2304 (1995). 



